Phonon-induced relaxation of a two-state system in solids 



O 
O 

(N 



Jaroslav Albert, E. M. Chudnovsky, and D. A. Garanin 

Physics Department, Lehman College, City University of New York, 
250 Bedford Park Boulevard West, Bronx, New York 10468-1589, U.S.A. 
(Dated: February 6, 2008) 

We study phonon-induced relaxation of quantum states of a particle (e.g., electron or proton) 
in a rigid double- well potential in a solid. Relaxation rate due to Raman two-phonon processes 
have been computed. We show that in a two-state limit, symmetry arguments allow one to express 
these rates in terms of independently measurable parameters. In general, the two-phonon processes 
dominate relaxation at higher temperature. Due to parity effect in a biased two-state system, their 
rate can be controlled by the bias. 

PACS numbers: 03.65.Yz, 66.35.+a, 73.21.Fg 



43 

o 

CD 



CO 



I 

C 

O 

o 



> 

\o 
\o 
in 
o 
\o 
o 



-a 
c 

o 
o 



X 



I. INTRODUCTION 

Relaxation and decoherence in a two state system cou- 
pled to environment is a fundamental problem of quan- 
tum physics. It was intensively studied in the past; see, 
e.g., the review of Leggett et al. lj. In the last years the 
interest to this problem has been revived by the effort 
to build solid state qubits. Recently, symmetry implica- 
tions have been considered for the problem of a particle 
in a rigid double- well potential embedded in a solid 0- 
It was demonstrated that symmetry arguments allow one 
to obtain parameter-free lower bound on the relaxation 
of quantum oscillations in a rigid double well, caused by 
the elastic environment. One of the arguments is that the 
double- well potential formed by the local arrangement of 
atoms in a solid is defined in the coordinate frame of that 
local atomic environment, not in the laboratory frame. 
Another argument is that interactions of the tunneling 
variable with phonons must be invariant with respect to 
global translations and rotations. When these arguments 
were taken into account, a simple universal result for the 
relaxation rate was obtained [2j in terms of measurable 
constants of the solid, with no unknown interaction con- 
stants. 

The above mentioned universal result refers to the low 
temperature limit when the relaxation of a two state sys- 
tem is dominated by the decay of the excited state due 
to the emission of one phonon. In this paper we ex- 
tend the method developed in Ref. 2 to the study of 
two-phonon Raman processes in double-well structures 

SHHISIili such processes can dominate relaxation 

at higher temperatures 0, IToL ITlj . We will show that 
in the temperature range bounded by the level splitting 
from below and by the Debye temperature from above, 
the rate of the Raman process for a biased rigid double 
well is given by a universal expression, very much like the 
rate of the direct one-phonon process. The Raman rate is 
proportional to the seventh power of temperature, while 
the one-phonon rate is linear on temperature. Interest- 
ingly, however, at small bias, the Raman rate, unlike the 
one-phonon rate, is proportional to the square of the bias. 
Consequently Raman processes can be switched on and 
off by controlling the bias. This universal result, which is 



a consequence of the parity of quantum states, must have 
important implications for solid-state qubits at elevated 
temperatures. Indeed, for an electron in a quantum dot, 
the rate of a direct one-phonon process is usually small. 
If the rate of a two-phonon process can be made small 
as well, this means that one can eliminate phonons as 
a significant source of relaxation and decoherence of the 
electron states in solid-state qubits. 



II. PARTICLE AND PHONONS 

Throughout this paper we shall use units where h = 
k,B = 1 unless stated otherwise. In the absence of 
phonons, the Hamiltonian in the laboratory frame is 



(1) 



where r is the radius vector, p is the momentum, and m 
is the mass of the particle (e.g., electron). A long- wave 
phonon described by the displacement field u(r) trans- 
lates the rigid double well in space. The Hamiltonian of 
the system (including the free phonon field) in the labo- 
ratory frame becomes 



H = f - + V(r - u) + H ph 
2m 



(2) 



Here, H. p h is the Hamiltonian of the free phonon 
field. We intend to obtain a Hamiltonian of the form 
H = H + Hph + He-ph where the last term describes 
the interaction of phonons with the electron in the dou- 
ble well potential. Using the fact that u is small, one can 
expand V(r — u) in Taylor series to obtain 



n, P 2 „/ S n. dV 1 9 V 



UiUj 



(3) 



The first three terms form the interaction-free part of 
the total Hamiltonian. The rest of the terms containing 
powers of u comprise the electron phonon interaction. It 
is clear that Eq. J2J requires detailed knowledge of the 
potential and its derivatives. One can, however, obtain 
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Eq. J5J) by performing a unitary transformation on Eq. 
JTJ with the help of translation operator TZ— e lp u 



n = e -ip.u noe i P u +Hph 

This can be expanded for small u as 
H = H + H p h + 

i 2 

+ i[H ,Pi]ui + —[[H ,Pi],Pj}ui 



(4) 



(5) 



Working out the commutators brings one back to Eq. 
©. However, the use of Eq. JSJ that we are going to 
employ allows one to obtain parameter-free results solely 
in terms of the energy levels of our effective two-state 
system without knowledge of the explicit form of V(r). 

We consider the case in which the particle, with good 
accuracy, is localized near r = ±R Q , where ±R D are the 
energy minima of the left or right wells. Without loss of 
generality we assume that ±R Q = ±X Q e x . The localiza- 
tion length of the state inside each well is small compared 
to the distance between the minima of the double-well 
potential. The bare ground states (when tunneling is ne- 
glected) in the left and right wells, that we denote by 
| ± X Q ), are approximately orthonormal, 



(±X \±X ) = l, (-X o \X o ) = 0. 



(6) 



The tunneling between the wells leads to the hybridiza- 
tion of the states given by orthonormal wave functions 



|0d 



V2 



(C±\X )tC t \-X )) 



where 



C± = y/l±e/A, 



VAT 



(7) 



(8) 



with A Q being the tunnel splitting in the unbiased double 
well and e being the energy bias between the wells. Note 
that the double well also has states, | with energies, 
Eg, other then E± corresponding to \ip±). The energy 
splitting 



A = E+ - E_ 



(9) 



is considered small compared to the distance from E± 
to other Eg. As we shall see, in this limit the summa- 
tion over all states |^) renders result for phonon-induced 
transitions between \ip±) that is insensitive to the explicit 
form of the potential. 

Below we shall deal with the matrix elements of op- 
erators p = p x , x, and their combinations. Other com- 
ponents of p and r are irrelevant. Localization of |^>±) 
allows one to compute matrix elements of powers of the 
operator x with the help of the relation 



x\±X ) =±X \±X ) 



(10) 



This gives 
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(V>±|a#±) = X ~(C 2 + - C 2 _) = X (e/A) 
(^-|a#+) - X C+C- = X (A /A) 
0>±\Ai>±) = *l> =0. (11) 

To compute other matrix elements we shall use relations 



and 



p = mx = —im[x^'H \ 



dx 2 



(12) 



px + xp — rnixx + xx) — m = —im[x 2 1 'H \ . (13) 

at 

This gives 

(^f = im(E e - E£,){ipe\x\ip?) 

(^\px + xp\ip e ) = im(E 6 - E e )(^\x 2 \ip 6 ,) (14) 

and thus, (ip-\px\'<p+) = 0. 

As we shall see, perturbation theory for Raman pro- 
cesses requires computation of the sum 



Ec — E-L 



(15) 



Application of Eq. I|12(l eliminates the denominator and 
yields 

E = im'J2{ip-\p\<ip i ){ipt\x\ip + ) . (16) 

Using the completeness of |^) we obtain 

S = im[(rl>- \ P x\iP+) - (V- \ P \^+) (V+ |ar|^+)] . (17) 

Finally, with the help of the above relations for matrix 
elements of x, p, and px we get 



E = -m 2 X 2 {A /A)e. 



(18) 



This is a mechanism of elimination of unspecified energy 
levels Eg from the problem, leading to a universal result. 



III. RAMAN MATRIX ELEMENT 

We are interested in the transition rate between the 
eigenstates of H Q + "Hph 



|V±) ® \<t>±) ■ 



(19) 



Here, \ip±) are the tunnel split states of the double well 
given by Eq. Q. The states \<j>±) are the eigenstates of 
TLph with energies E p h.±- Our goal is to study Raman 
scattering processes involving absorption of a phonon 
of frequency Wk and emission of a phonon of frequency 
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FIG. 1: Raman process on tunnel-split levels in a double- well 
potential, including virtual transitions to higher levels [the 
first term in Eq. 1241 shown]. 



= + A, accompanied by the transition of the par- 
ticle \ip+) — > \t(>-) , see Fig.l. The Raman rate can be 
computed with the help of the Fermi golden rule in the 
second order in the interaction. The matrix element for 
this process is the sum of two matrix elements, 

M = M {2) +M (1+1) . (20) 

The first term denotes the first order perturbation on 



H 



(2) 

e—ph 



i 



[[H ,p],p]ul, 



(21) 



while the second term stands for the second order per- 
turbation on 



^e-ph = i[H-o,P]u x . 

The explicit expressions for and M^ 1+1 ^ are ^ 



(22) 



and 

M (i+D 



M (2 



E 
E 



(2) 



-ph l*+ 



(23) 



(i) 



E+ + fitJk — E$ 



E + — HlUq — 



(1) 



(24) 



Here, \^>^) is a direct product of the eigenstates of H 
with the phonon states |nk — 1, fi q ) in the first term and 



l^kj + 1) m the second term. 



First, we calculate the phonon parts of and 
m( 1+1 ) using canonical quantization of the phonon dis- 
placements 



u 



1 ^e kA e lkr 



— T 



/^kA 



(25) 



where p is the density of the solid, V is its volume, ekA 
are unit polarization vectors, A = t\,t2,l denotes polar- 
izations, and WkA = v\k is the phonon frequency that we 
will usually write as Wk- Writing \^f+) = |rak,rc q ) (g> 



and l^-) = |nk — 1, n q + 1) ® |^_) , for the phonon ma- 
trix element in one obtains 

(n k - l,n q + 1| u 2 |n. k! n q ) 



= 4 e ^ (k - q) - r W nkK+1) ^M pfe . (26) 

For the phonon matrix elements in M one obtains 

(n k - l,n q + 1| |n k - l,n q ) (n k - l,n q | |n k ,n q ) 
= (nk - l,n q + 1| |n k ,n q + 1) (n k ,n q + 1| u x |n k ,n q ) 
= M ph /2, (27) 

(One can see from the completeness relation that the sum 
of these two expressions is M p h ) 

Next, we evaluate the particle parts of and 
A/U+i). For (V>_| [[H p],p\\tp + ) that enters M^ 2 \ writ- 
ing the commutator explicitly and inserting the identity 
operator 1= |^) (^| results in 

(1>- 1 (Hop 2 - 2pH oP + p 2 H ) |^+) 
= ^(£+ + E_ - 2E t ) (</;_ | p .(28) 

e 

The particle part of M' 1+1 ) simpifies to 

= (E i - -£+) (Ve|P^+) -(29) 

Collecting the terms, for M of Eq. (g^l one obtains 



(30) 



where 



Q i = E+ + E_- 2E S + {E 6 - E-){E t - E+) 
1 1 



E, 



(31) 



A common mistake that propagates through literature 
[l4j is summation of rates due to M^ 2 " 1 and Af( 1+1 ', in- 
stead of adding matrix elements first and then squaring 
the result and computing the rate. This mistake is not in- 
nocent since 

M<-V and M^ 1+1 ) may cancel leading parts 
of each other. Taking into account conservation of en- 
ergy w q = cjk + A and the relation = E + — A, one 
can rewrite this expression as 



2oj k K + A)(£ g - E+ + A/2) 
(£7 £ - - w k )(£? - ^+ + w k + A) ' 



(32) 



One has Q± = =pA. We consider the case A <C w k ~ 
T^; where u>d is the Debye frequency. Then for £ 7^ 
± the energy differences E^ — E+ are considered to be 
large compared to A, that yields 



= = gtJkjaj. + A) 



(33) 
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From Eq. fTty follows that the terms with £ = ± in Eq. 
(13011 disappear. Using Eq. (|15fl . one obtains 

m = \m p}1 y, y>-\p\4>t)y>t\p\>i>+)Q<i 

= M p/l w k K + A)E (34) 
and finally, with the help of Eqs. (ffi>|) and ^18(1. 

„ 1 x x /"k(n q + l) 2v2 A Q 

M = u e kA e qeA/ — ^ J -m 2 X 2 e—uj k w k + A . 

pV 4 y WkWq A 

Here we have suppressed an irrelevant phase factor. This 
result for the Raman matrix element is insensitive to the 
explicit form of the double- well potential V{x). It would 
be a hopeless task to obtain it from Eq. (J3J. 

IV. RAMAN TRANSITION RATE 



According to theory of elasticity 0] vi > \[2v t . Thus, 
the second term in this expression is a small correction 
and it can be neglected. Since we are interested in the re- 
gion A <^T we can keep only the leading order on T/A in 
Eq. (|3"9"|) . The Raman rate for A <CT^C uj p> then becomes 

where we have introduced characteristic energy and fre- 
quency scales 

£ «= m " ^Gs&r <42) 

of the problem that are entirely determined by the pa- 
rameters of the unperturbed dot and its elastic environ- 
ment. 



According to the Fermi golden rule |l2j the Raman 
rate is given by 



E 



V 2 \M\ 2 2ir5(u>q — cjk — A). 



(36) 



The integration variables dkdq can be written in spheri- 
cal coordinates as dkdq k 2 q 2 dfl^dH.q, and the integral for 
T<ti lod can be expressed as 



r 2 = 



where 



2w 3 p 2 

x (a) 



o 



B = 



Att 2^1 



dLUkdujqUJkUJqn^nq 
5(LJq - CJ k - A), 



-1) 

(37) 
(38) 



v\ is the velocity of sound with polarization A, whereas 



n k 



l) is the Bose occupation number of a 



phonon. For Wk ~ T 3> A one obtains the Raman rate 



m 4 X Q 4 A 5 Agg z 
2^p 2 h 11 



1 + 



(39) 



where a = J °° dxx 6 e x (e x — 1) 2 = (16/21) it 6 . One can 
see from this integral that the validity of Eq. (|39|l re- 
quires at least T < cjd/10. B can be calculated with the 
help of the transverse-phonons sum rule Y]f— t t (ekt ■ 

a)(eki • b) = (a ■ b) — (a ■ k)(k • b), where k = k/fc, etc. 
Setting a = b = and averaging over the directions of 
k yields 



V. DOT FRAME CALCULATION 

In this section we will check our result by calculating 
the Raman rate in the frame of reference of the dot, as it 
was done for one-phonon processes Q . In the laboratory 
frame the Lagrangian of the particle is 



r m c' 



u 



- V(r') 



(43) 



where r' is the radius vector of the particle of mass m in 
the coordinate frame rigidly coupled to the double well. 
The linear momentum that is canonically conjugated to 
r' is given by 

Br 

P' = ^T = m (*' + *) ■ ( 44 ) 
or 

The corresponding Hamiltonian is 

H P (p', r') = p' • r' - C P = ^- - p' • u + U(r') . (45) 

The full Hamiltonian is Tip + TL p h ■ Contrary to the pre- 
vious model described by Eq. ©, we now have only one 
interaction term, — p' • u. Similarly to Section III, one 
can write the matrix element for the Raman processes as 



M 



E 
E 



(-0-1 p ■ u|-0c) (AIp • "IV'H 

E+ + fio^k — 
(V'-Ip' -u|V>£) (^Ip' •uIV'h 



E-L. — frhJct 



E : 



(46) 



Inserting [lj; 



B 



f dQ k 


(I" 




1 4tt 







1 



(40) 



— E 

2pV ^ 



l -kA , 



kA 



(47) 
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into Eq. (|46|l and evaluating the matrix elements, we ob- 
tain 

s 

which coincides with Eq. (|30|) up to an insignificant 
phase. 

VI. DISCUSSION 

We have demonstrated that the two-phonon relaxation 
of the tunnel-split states of a particle in a biased solid- 
state double-well potential can be expressed in terms of 
independently measured parameters and without any un- 
known constants. Two-phonon processes may dominate 
relaxation at elevated temperatures (see below). An in- 
teresting observation, however, is that at a small bias the 
rate of Eq. Q41|l is proportional to e 2 , while at a large 
bias it becomes independent of e. This means that one 
can switch Raman processes on and off by controlling the 
bias. This result may seem strange at first, however, it is 
a fundamental consequence of quantum mechanics. The 
reason for this effect is parity. If we remove the bias, 
the potential well will become symmetric. Consequently, 
the Hamiltonian and the parity operator commute which 
leads to eigenstates of even or odd parity. It is easy to see 
that the states and at e = have even and 

odd parity, respectively. Therefore, the matrix elements 
in Eq. will all vanish. 

To find the range of parameters where two-phonon re- 
laxation becomes important, the rate of the Raman pro- 
cesses should be compared with onc-phonon transition 
rate Q that can be written in the form 

Notice that Eqs. (|49l) and (|49|l do not contain any un- 
known interaction parameters. The quantity of interest 
is the ratio T2/T1 which can tell us the importance of the 



second order process versus the first order process at var- 
ious temperatures. If we take T > A, the coth(A/2T) in 
Eq. (|41[) can be replaced bv 2T/A. The above mentioned 
ratio then yields 




At any given temperature this ratio has a maximum at 
e = A . For, an electron in a double- well dot with 
X Q ~ 10 nm embedded in (or deposited on) a solid with 
P ~ 5 g/cm 3 and vt ~ 10 3 m/s, parameter £ t is of order 
300 K. Then, for, e.g., £ ~ A„ ~ 1 K, Raman processes, 
according to Eq. (|50|l , will dominate electron-phonon re- 
laxation above 33 K, while below that temperature the 
relaxation will be dominated by direct processes. The 
actual phonon rates for an electron are not likely to ex- 
ceed 10 6 s _1 even at T ~ 100 K. For a proton in a 
molecular double well with X a ~ 0.3 nm in a solid with 
p ~ 5 g/cm 3 and Vt ~ 10 3 m/s, one gets Et ~ 40 K. At 
£ ~ A Q ~ 1 mK, according to Eq. JSUJ), Raman pro- 
cesses will dominate proton-phonon relaxation above 1 
K, while direct processes will dominate relaxation in the 
millikelvin range. 

Finally we should note that since our model is based 
upon bare quantum states that are well localized in space, 
it is rigorous for heavy particles, like, e.g., a proton or 
an interstitial atom, but is less rigorous for such a light 
particle as an electron. Nevertheless, even for an elec- 
tron our formulas should provide a good approximation 
in the limit of weak tunneling between the wells. Note 
also that at a large tunnel splitting, the actual rates for 
a heavy particle like proton, interstitial atom or defect, 
can become so large that the approximation based upon 
Fermi golden rule will no longer apply |16| . Even in this 
case, however, the matrix elements can be expressed in 
terms of measurable parameters of the quantum well and 
the solid. 
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